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Cooperative Multiscale Aging in a Ferromagnet/Antiferromagnet Bilayer 
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We utilize anisotropic magnetoresistance to study temporal evolution of the magnetization state in 
epitaxial NisoFe 2 o/Fe 5 oMn 5 o ferromagnet/antiferromagnet bilayers. The resistance exhibits power- 
law evolution over a wide range of temperatures and magnetic fields, indicating that aging is char¬ 
acterized by a wide range of activation time scales. We show that aging is a cooperative process, 
i.e. the magnetic system is not a superposition of weakly interacting subsystems characterized by 
simple Arrhenius activation. The observed effects are reminiscent of avalanches in granular materi¬ 
als, providing a conceptual link to a broad class of critical phenomena in other complex condensed 
matter systems. 

PACS numbers: 85.70.Kh,89.75.-k,89.75.Da 


In bilayers of materials with different lattice parame¬ 
ters, structural frustration can result in dislocations or 
even amorphous interlayers [ij] . Similarly, magnetic frus¬ 
tration can be expected at interfaces between materials 
with different magnetic orders Q. The origin of the frus¬ 
tration is the random effective field experienced by both 
materials due to their exchange interaction across the in¬ 
terface that is generally not atomically smooth. In partic¬ 
ular, some of the unusual magnetic properties exhibited 
by bilayers of antiferromagnets (AF) and ferromagnets 
(F) have been attributed to the magnetic domain walls 
that are formed to reduce the interfacial exchange en¬ 
ergy S3 , or even disordered spin states near the F/AF 
interface @, @|. After almost 60 years of extensive re¬ 
search, fundamental understanding of F/AF bilayer sys¬ 
tems remains elusive. Besides theoretical challenges in 
describing the effective exchange fields at the F/AF in¬ 
terfaces, common experimental approaches, such as the 
hysteresis loop measurement, can lead to irreversible 
changes of the magnetic configuration, thus obscuring 
the essential signatures of frustrated systems such as ag¬ 
ing m. Therefore, nonperturbative real-time charac¬ 
terization methods may be needed to provide insight into 
the properties of these systems. 

We utilized anisotropic magnetoresistance (AMR) to 
characterize the evolution of the magnetization state in 
F/AF bilayers without perturbing the system. Our mea¬ 
surements reveal power-law relaxation over a wide range 
of temperature T , indicating multiple scales of activa¬ 
tion energies, and activation times spanning an estimated 
range of at least seven orders of magnitude. These results 
are not affected by the variations of field H, demonstrat¬ 
ing that activation occurs in the AF layer. The observed 
power-law form of relaxation is also independent of the 
prior magnetic aging history, indicating that aging is a 
cooperative process; it cannot be described in terms of 
independent activation barriers. These results provide 
an unprecedented insight into the magnetism in F/AF 
bilayers, linking them to other complex and frustrated 
systems that exhibit cooperative aging phenomena. 

Our samples were deposited by high-vacuum sput¬ 


tering on (OOOl)-oriented sapphire substrates annealed 
in air at 1300° C to achieve atomically flat surface. 
An epitaxial (lll)-oriented Pt(5) buffer layer was de¬ 
posited at 550° jsj, followed by the F/AF bilayer 
Ni 8 oFe 2 o( 10 )/Fe 5 oMn 5O ((i) deposited at room tempera¬ 
ture to avoid interdiffusion of the magnetic interfaces. 
All thicknesses are in nanometers (nm). The bilayer was 
capped with SiO2(20) to prevent oxidation. We have fab¬ 
ricated and studied several samples with thickness d of 
FeMn ranging from 1 to 3.5 nm. Their magnetic prop¬ 
erties, such as the temperature Tb characterizing the 
onset of h yst eresis loop asymmetry - the exchange bias 
(EB) flQ, ll|, were consistent with the previous stud¬ 
ies jl2j - ll5l ] . We focus on a sample with d = 2 nm, which 
exhibited Tb = 140 K within the range of temperatures 
T = 5 — 300 K accessible in our magnetoelectronic mea¬ 
surements. We will show below that some of the signa¬ 
tures of EB extend far above Tb . 

The magnetoelectronic characterization was performed 
in the four-probe van der Pauw geometry, using an ac 
current / = 0.1 rnA rms at frequency / = 1.3 kHz, 
and lock-in detection of ac voltage V. The resistance 
R = V/I exhibits a 180°-periodic sinusoidal dependence 
on the in-plane orientation of a sufficiently large in-plane 
field H [Fig. [lja)], as expected due to the AMR of the 
Permalloy(Py)=Ni 8 oFe 2 o layer. The angle 0 = 0 corre¬ 
sponds to the held orientation perpendicular to the cur¬ 
rent. The monotonic increase of the magnetoresistance 
A R = i?(90°) — -R(O) with decreasing temperature T is 
not affected by EB [Fig. [2(b)], confirming that the mag¬ 
netoresistance is determined entirely by the AMR of Py. 

When H is swept at 6 = 0, the resistance exhibits 
sharp peaks at the coercive fields H + , H_ [Fig.QJc]]. To 
establish EB, we used the conventional procedure [lOj, 11] 
of cooling from room temperature RT = 295 K in satu¬ 
rating held H = 500 Oe, at 0 = 0. The hysteresis loop 
becomes asymmetric below Tb = 140 K, as illustrated 
in Fig. [2(c) for T = 5 K. We note that the magnetore¬ 
sistance (MR) observed in the hysteresis loop is close to 
A R determined from the rotational AMR [Fig. [2(b)], in¬ 
dicating that reversal occurs through the conhguration 
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FIG. 1: (Color online) Magnetoelectronic characterization of 
the F/AF bilayer: (a) Dependence of resistance on the in¬ 
plane orientation of the field, at T = 295K, H = 100 Oe (solid 
symbols), and 5 K, H = 2 kOe (open symbols), (b) Tempera¬ 
ture dependence of magnetoresistance determined from R(9) 
data obtained at H = 2 kOe, (c) Symbols: magnetoelectronic 
hysteresis loop obtained at 5 K, after 4 prior similar ’’train¬ 
ing” loops, up(down) triangles are for increasing(decreasing) 
H. Curves: fits based on the effective exchange field model, as 
described in the text. Coercive fields and H+ are labeled, 
(d) Temperature dependence of coercivity He = (H+ — H -)/2 
and exchange bias field He = ~{H- + H+)/2. 


FIG. 2: (Color online) Magnetic aging at T = 5 K. (a) Sym¬ 
bols: time evolution of resistance R at Hf = —300 Oe, after 
prior aging at Hi = 350 Oe for 1000 s. Curve: fitting with the 
power-law dependence R(t) = Ro + At~ c , with Ro = 3.496 D, 
A = 0.011, c = 0.081. Inset: R—Rq vs t plotted on the log-log 
scale, (b) Seven sequential aging cycles at Hf = —500 Oe, 
as labeled, each preceded by aging at Hi = 350 Oe for 1 Ks. 
(c,d) Dependence of the exponent c (c) and the asymptotic 
resistance Ro (d) on the aging cycle. The measurements were 
performed immediately after cooling at H = 500 Oe from 
RT. The shown fitting error bars underestimate the uncer¬ 
tainty of the determined values, due to the additional errors 
from transients at small t. 


of M almost homogeneously transverse to the field, con¬ 
sistent with a high spacial uniformity of the magnetic 
properties. The dependencies of both the EB field He = 
— ( Hi + H 2 )/2 and the coercivity He = (H 2 — H\)/2 on 
T are consistent with other studies of EB systems based 
on FeMn [lj-[ijl- 

Our central result is the demonstration of cooperative 
multiscale aging over a wide range of T and H. To ob¬ 
serve aging, the field H was ramped at a rate of 2 kOe/s 
from the initial value Hi above H + [or below if_] to 
a final value Hf below H_ [or above H+], and subse¬ 
quently R was recorded in 1 s time increments, as illus¬ 
trated in Fig. d a ) for T = 5 K, Hf = —300 Oe. The 
lock-in time constant was set to 100 ms to minimize the 
effects of the instrumental bandwidth. We emphasize 
that both H and T were constant during aging, i.e. the 
recorded evolution was not perturbed by the measure¬ 
ment. Thermal activation in F/AF heterostructures has 
been extensively discussed in the context of granular sys¬ 
tems, where the particle size is expected to set the en¬ 
ergy scale Eg for the Arrhenius-type exponential decay 
R(t) = Rg + i?i exp[—f/r[ over the characteristic time 
r oc exp[Eg/kT] [16j, LL7| - Here, k is the Boltzmann con¬ 
stant. In our measurements, the evolution of R(t) for t 


up to about 50 s could be well fitted with simple expo¬ 
nential decay. However, such fitting became inadequate 
at longer time scales. To remedy this discrepancy, one 
can assume a certain distribution of activation barriers, 
providing additional fitting parameters [l8|. In the ex¬ 
treme limit, there is no characteristic activation energy 
scale, and consequently R(t) does not exhibit a charac¬ 
teristic decay time. A similar situation is encountered at 
critical points in phase transitions, resulting in power-law 
dependencies of physical properties [19j. Indeed, power- 
law dependence R = Ro + At~ c provided an excellent fit 
for all of our experimental R(t) data, as illustrated by the 
curve in Fig. da) for data spanning three orders of mag¬ 
nitude in t. The power-law dependence was observed at 
temperatures from 5 K to over 200 K, indicating that the 
activation energies are spread over more than two orders 
of magnitude. 

The aging curves depended on the previous aging his¬ 
tory, reminiscent of the training effect - variations of the 
magnetic properties observed in sequential magnetic hys¬ 
teresis loops of F/AF bilayers 5. |20|. Figure db) shows 
seven aging curves acquired at Hf = —500 Oe. Each 
curve was measured after pre-aging at Hi = 300 Oe over 
time At, = 1 Ks. The overall form of the dependence 
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FIG. 3: (Color online) (a) R vs t at T = 5 K, Hf = 150 Oe, 
after prior aging at Hi = —300 Oe over time interval At; = 
500 s (bottom curve), 200 s, 100 s, 51 s, 21 s, 11 s, 6 s, and 3 s 
(top curve). (b,c) Power-law exponent c (b) and relaxation 
scale A (c) vs A ti (symbols), and a logarithmic fit to the data 
(line). Note the logarithmic scale for A ti. 


FIG. 4: (Color online) Dependence of the aging character¬ 
istics on field, at T = 5 K. (a) The power-law exponent c 
determined from fitting the Rvst. data. Error bars show the 
fitting uncertainty, (b) Symbols: experimental dependence 
of the total relaxation amplitude R(t = Is) — Ro on field. 
Lines: fitting based on the effective exchange field model, as 
discussed in the text, using the data for \H\ > 200 Oe. 


remains similar for different aging cycles, as confirmed 
by fitting with the power-law dependence, with expo¬ 
nent c that exhibits only a modest irregular dependence 
on the aging cycle [Fig. 0(c)] . The main difference be¬ 
tween the consecutive aging curves is the overall increase 
of R. Since the value of R itself evolves in time, the 
increase of R in consecutive cycles can be characterized 
by its asymptotic value Rq , which exhibits a monotonic 
increase with the aging cycle number [Fig. [2][c)] . This re¬ 
sult demonstrates that the state asymptotically reached 
by the magnetic system after a single aging cycle is de¬ 
pendent on the magnetic history, similar to the kinetic 
trapping in glasses [21, 22], The data of Fig.[2]confirm the 
dependence of the magnetic properties, including aging 
characteristics, on the specific protocol used to prepare 
the magnetic system, which is well known both in the EB 
community and in studies of glassy systems. The results 
described below were obtained after multiple aging cy¬ 
cles, in the regime where the main aging characteristics 
stabilize [see Figs. [2jc,d)]. We anticipate that studies of 
the dependence on the preparation protocol may provide 
insight into the aging mechanisms. 

The dependence of i?o on the aging history indicates 
that the system cannot be described in terms of indepen¬ 
dent activation barriers. To further test this conclusion, 
we performed aging measurements in which we varied 
only the time At r of pre-aging in the reversed state, from 
500 s down to 3 s [Fig. 0(a)] • The shape of the aging 
curves did not depend on At,;, and the decay exponent c 
obtained from the power-law fitting remained the same 
within the approximately 10% data spread [Fig. [3jb)] . 
If the relaxation could be described by independent ac¬ 
tivation processes characterized by the barriers E n and 
the corresponding relaxation times t„, 1 < n < N, then 
the subsystems with Tn ^ A ti would become activated, 
while the subsystems with r n > At, would not be acti¬ 
vated. As a consequence, for small A ti the aging curves 
would exhibit significantly smaller amplitude of decay at 
t > A ti. Since the form of the experimental aging curves 


is independent of A ti, we conclude that aging involves co¬ 
operative processes coupling multiple energy scales, and 
cannot be described by the Arrhenius-type activation of 
individual weakly coupled subsystems. 

While the exponent c characterizing the aging curves 
remained independent of Ati, the decay scale A exhib¬ 
ited a small but well-defined decrease by about 15% 
when Ati was decreased by two orders of magnitude 
[Fig. 0(c)]. The values of Ati below a few seconds are 
not experimentally accessible in our magnetoelectronic 
measurements. Nevertheless, one can generally expect 
that A should vanish when Ati becomes smaller than 
the shortest activation timescale in the system. The ap¬ 
proximately logarithmic dependence in Fig. He) extrap¬ 
olates to an intercept A = 0 at 10 -4 s. Thus, activation 
timescales likely span at least seven orders of magnitude, 
from 10~ 4 s to at least our characteristic measurement 
time of 10 3 s. We note that the extrapolation of aging 
characteristics to t = 10” 4 s does not lead to unphysi¬ 
cal results. For instance, the total estimated variation of 
resistance, R(t = 10 _4 s) — Rq = 10 4c A, did not exceed 
0.12 fl in all of the aging experiments at 5 K. This value 
is smaller than the full MR A R = 0.28 12 at 5 K. 

Both F and AF can contribute to aging in F/AF bilay¬ 
ers. The magnetic anisotropy and/or defects in F and/or 
AF can define the energy barriers for the magnetic evo¬ 
lution, which can proceed either by the uniform reversal 
in small magnetic grains, or by F or AF domain wall mo¬ 
tion in continuous films. To establish the relative con¬ 
tributions of the two magnetic layers, we determined the 
dependence of aging characteristics on the field Hf. The 
Zeeman energy contribution should result in the expo¬ 
nential dependence of activation in the F layer on Hf, 
while the corresponding dependence for AF should be 
weak. Figure Ha) summarizes the values of the expo¬ 
nent c determined from aging at different Hf ranging 
from —650 Oe to 650 Oe. The magnetoelectronic signa¬ 
tures of aging became too small for reliable measurements 
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at | Hf > 650| Oe [see Fig.0[b)]. The values of c exhibit 
random variations around the average c = 0.084, and no 
correlation with Hf. Based on this result, we conclude 
that aging involves activation processes in AF that af¬ 
fect the magnetization M of Py only through exchange 
coupling at the interface and are not directly influenced 
by H. The error bars in Fig. 0a) reflect only the fitting 
uncertainty, and do not account for the additional er¬ 
rors caused by the transient effects at short times scales, 
which are caused by the limited bandwidth of both the 
electromagnet power supply and the lockin amplifier. An 
additional measurement time error is caused by the onset 
of aging while the field is still being ramped. 

While the form of the aging curves was independent 
of Hf, the overall scale of the resistance decay rapidly 
decreased with increasing magnitude of Hf, as shown by 
symbols in Fig0b). Here, we plot the difference between 
the first measured resistance value at R(t = Is) and its 
asymptotic value Rq , which, as we shall see below, is more 
convenient than the scale A for the quantitative analysis 
of relaxation. The dependence in Fig. 0b) is consistent 
with our conclusion that aging occurs in the AF layer. 
We can describe the exchange interaction of F with AF 
by an effective exchange field H' with average compo¬ 
nents ifj|, H' L in the direction of H and perpendicular 
to it, respectively. Both of these components vary over 
time due to the AF aging, resulting in variations of the 
angle <f> formed by the Py magnetization M relative to H, 
according to (/> « H' ± /(H + i?|), and the corresponding 
variations of resistance 


R — Rmin T 2 

where Rmin is the resistance minimum at </> = 0. Our 
model is supported by the excellent agreement of the 
fit based on Eq. 0 with the quasi-static measurements 
of R vs H [solid curves in Fig. [Q[c)], yielding H' ± = 
148.5 ± .5 Oe from two independent fits of both hystere¬ 
sis loop branches up to the switching points. Fitting the 
same curves beyond the switching point is less meaning¬ 
ful, because of the aging that occurs concurrently with 
the field sweep. Using the form of Eq. © to analyze 
the dependence of the relaxation magnitude on Hf in 
Fig. 0b), we obtain a good fit for all the \H\ > 200 Oe 
data, with a single set of fitting parameters H'^ = 50 Oe, 
A H' ± = 105 Oe [curves in Fig. 0b)]. Here, A H' ± is the 
overall reduction of the effective transverse field due to 
aging between t = Is and oo. Thus, the dependence 
of the aging curves on Hf can be explained entirely by 
the effect of H on M, while its direct effect on aging is 
negligible. 

Behaviors similar to those discussed above for T = 5K 
were observed at higher T, extending significantly above 
the blocking temperature Tb = 140 K [see Fig. [Sja) for 
aging curves and their power-law fits at 210 K]. To quan- 


H’ ± 


H + m 


(i) 





FIG. 5: (Color online) Dependence of aging characteristics on 
temperature, (a) Symbols: measured R vs t at T = 210 K, 
Hf = —65 Oe and Hf = 55 Oe, as labeled. Curves: fits 
with the power-law dependence, and power-law exponents c = 
0.126 and 0.246. (b) Hf used in the measurements of aging 
at different T. (c,d) Power-law exponent c (c) and scale A (d) 
vs T, obtained from fitting of the aging data such shown as 
in panel (b). Filled (open) symbols are for Hf > 0 (Hf < 0). 


titatively characterize the dependence on T, aging was 
measured at the values of Hf adjusted so that the cor¬ 
responding resistance in the hysteresis loop was approx¬ 
imately at 20% of A R above the minimum [Fig. 0(b)]. 
The power-law exponent c exhibited an overall increase 
from 0.08 at 5 K to about 0.2 around 200 K [Fig. 0(b)] . 
This trend is superimposed with increased random vari¬ 
ations of c, which can be correlated with the decrease of 
the relaxation scale A above Tb [Fig. 0(d)], preventing 
reliable measurements of aging at T > 210 K. The in¬ 
crease of c with increasing T is consistent with the larger 
relative contribution of fast activation processes, and a 
smaller contribution of slow processes. Nevertheless, the 
power-law form of aging curves is retained up to 210 K 
[Fig. 0a)[, indicating non-negligible contribution from 
the activation times r ss 1000 s. The decrease of A at 
high T can be attributed to the decreasing volume of 
AF involved in aging on the experimentally accessible 
timescales. 

Summarizing our main findings, magnetoelectronic 
measurements reveal aging in epitaxial thin-film 
F/AF=Py/FeMn bilayers over a wide range of tempera¬ 
tures. The aging curves exhibit simple scaling with the 
applied magnetic held, demonstrating that the activa¬ 
tion processes are confined to the AF layer. Aging is 
characterized by the power-law dependence on time and 
a weak dependence on temperature, indicating a wide 
range of the activation energies, and a range of the acti¬ 
vation times that at T = 5 K is estimated to extend over 
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at least seven orders of magnitude, from 10 -4 s to at 
least the longest measurement time of 10 3 s. Finally, the 
form of the aging curves is independent of the magnetic 
pre-aging history, indicating that activation processes are 
cooperative, in other words the magnetic system cannot 
be described as a superposition of weakly-coupled sub¬ 
systems activated according to the Arrhenius law. 

The observed power-law evolution is reminiscent of the 
self-organized criticality and the associated avalanche dy¬ 
namics [Hi- I n ferromagnets, such avalanches of mag¬ 
netic domain reversal are observed as the Barkhausen 


noise 24]. The emerging physical picture for aging in 
AF is that of stable AF regions interfaced with regions 
that can be either stable or activated, depending on the 
stability of their AF environment. Activation of a cer¬ 
tain region can activate or deactivate other neighboring 
AF regions due to their exchange interaction, which can 
result in avalanches of AF activation. According to the 
picture of self-organized criticality [HI, at long times the 
AF is expected to asymptotically form ’’minimially stable 
clusters”, which can be irreversibly perturbed by small 
variations of temperature or changes of the magnetic con¬ 
figuration of F. This picture implies a correlation between 
the effects of the magnetic and the thermal history. Test¬ 
ing this correlation can elucidate the mechanisms of co- 
operativity. 

The geometry of the activated AF clusters (and by ex¬ 
tension the stable AF regions) is likely fractal, since aging 
phenomena are driven by the random effective exchange 
held. Our hypothesis is supported by the apparent lack 
of the characteristic activation energy scale, since the lat¬ 
ter is determined by the AF anisotropy and scales with 
the volume of the activated regions. Further spatially re¬ 
solved studies, e.g. based on x-ray dichroism microscopy 
of AF and/or F layers will likely elucidate the spatial 
characteristics of activation [see Ref. [25] for microscopic 
images quite conducive of this physical picture]. Our re¬ 
sults may have significant implications for other F/AF 
heterostructures based on single-crystal and even poly- 
crystalline AF materials, where crystallinity may provide 
a natural limit for the geometry of the activated AF clus¬ 
ters and cooperativity, but some signatures of the behav¬ 
iors described above are likely retained. F/AF bilayers 
may also represent a controllable (by means of the mag¬ 
netic held and temperature) model system that can pro¬ 
vide insight into other condensed matter systems exhibit¬ 
ing complex critical phenomena. For instance, the ran¬ 
dom effective exchange held experienced by the AF can 
provide a straightforward implementation for the classic 
Imry-Ma problem of random-held magnetism [26|, [27]. 
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